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Abstract 

The  carbon  oxidation  rates  of  Nagle  and  Strick- 
land-Constable  are  routinely  applied  in  the  flame 
community  to  determine  soot  oxidation  by  molecu¬ 
lar  oxygen.  These  rates  were  experimentally  deter¬ 
mined  at  temperatures  from  1000  to  2000®C  (1273 
to  2273  K)  and  at  oxygen  partial  pressures  from 
0.1  to  0.6  atm  and  then  curve  fit.  For  typical  labora¬ 
tory  flames,  oxygen  partial  pressures  may  be  sig¬ 
nificantly  lower  than  0.1  atm.  Frequently,  research¬ 
ers  will  extrapolate  the  Nagle,  Strickland-Constable 
(hereafter  referred  to  as  NSC)  cunre  fits  to  these 
lower  oxygen  partial  pressures  to  determine  the 
molecular  oxygen  contribution  to  the  overall  soot 
oxidation  rate.  It  is  shown  that  the  NSC  curve  fits 
exhibit  nonphysical  behavior  at  pressures  below 
0.05  atm  (depending  on  the  temperature).  This 
behavior  is  revealed  when  the  oxidation  rates  are 
transformed  to  collision  efficiencies.  At  a  fixed  tem¬ 
perature,  the  oxidation  rate  drops  faster  than  the 
collision  rate  as  the  pressure  is  reduced.  At  typical 
flame  conditions,  extrapolation  of  the  NSC  curve 
fits  yields  collision  efficiencies  that  may  be  an  order 
of  magnitude  too  low. 

Introduction 

The  carbon  oxidation  rates  of  Nagle  and  Strick- 
land-Constable^  are  routinely  applied  in  the  flame 
community  to  determine  soot  oxidation  by  molecu¬ 
lar  oxygen.  These  rates  were  experimentally  deter¬ 
mined  at  temperatures  from  1000  to  2000°C  (1273 
to  2273  K)  and  at  oxygen  partial  pressures  from 
0.1  to  0.6  atm  and  then  curve  fit.  For  typical  labora¬ 
tory  flames,  oxygen  partial  pressures  may  be  sig¬ 
nificantly  lower  than  0.1  atm.  Frequently,  research¬ 


ers  will  extrapolate  the  Nagle  and  Strickland-Con¬ 
stable  (hereafter  referred  to  as  NSC)  curve  fit  to 
these  lower  oxygen  partial  pressures  to  determine 
the  molecular  oxygen  contribution  to  the  overall 
soot  oxidation  rate.  For  example,  Lin,  Sunderland, 
and  Faeth^  used  the  NSC  rate  expression  to  evalu¬ 
ate  soot  oxidation  by  molecular  oxygen  in  19-kPa 
(0.188  atm)  acetylene/air  diffusion  flames.  In  their 
flames,  the  oxygen  partial  pressure  ranged  from 
0.0004  to  0.004  atm  on  the  flame  axis.  Puri,  San¬ 
toro,  and  Smyth®  also  use  the  NSC  expression  to 
evaluate  soot  oxidation  by  O2  in  a  variety  of  flames 
at  oxygen  partial  pressures  in  the  neighborhood  of 
0.02  atm. 


Soot  Oxidation  Rates  and  Coiiision  Efficiencies 


Nagle  and  Strickland-Constable  fit  their  oxida¬ 
tion  rate  data  with  the  expression 


to  =12. 


V  +•^2%) 


ix  +  (l<BP02)(1-X)lgm  cm  sec 


;)jgm  cm  *sec 


where 


=  20exp{-15100/T),  gm  cm  ^sec  ^  atm  ^ , 
kg  =  0.00446exp(-7640/T),  gm  cm“^sec"^atm"^ , 

c  _ 2  — 1 

kT=  1.51  X  10  exp(-48800/T)gm  cm  sec  ,and 
kz  =  2.13exp(2060/T),  atm"’ . 


(1) 


The  specific  soot  oxidation  rate  (®)  is  given  as  a 
function  of  temperature  (T)  in  degrees  Kelvin  and 
oxygen  partial  pressure  (Pq  )  in  atmospheres. 
The  NSC  rate  expression  is  pitted  as  a  function  of 
oxygen  partial  pressure  and  lAemperature  in  Fig. 
1 .  The  expression  is  evaluated  only  over  the  exper- 
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imental  range  of  these  parameters,  i.e.  1273  to 
2273  K  and  0.1  to  0.6  atm.  Note  that  the  oxidation 
rate  is  a  monotonic  increasing  function  of  oxygen 
partial  pressure. 


Fig.  1 .  NSC  curve  fit  of  oxidation  rates. 

The  oxidation  rate  can  also  be  expressed  as  a 
fraction  of  the  totai  number  of  collisions  with  O2: 
i.e. 

CD  =  Zfm  (2) 


or 


Zifi 


(3) 


where  Z  is  the  coiiision  frequency,  Zf  is  the  fre¬ 
quency  of  reactive  collisions,  f  is  termed  the  coiii¬ 
sion  efficiency,  and  m  is  the  mass  removed  per 
reactive  collision.  The  collision  rate  with  O2  is  given 
by 

1  1  P.»  /sFT  , 

^  (4) 

where  k  is  the  Boltzmann  constant,  is  the  aver¬ 
age  molecular  velocity,  n  is  the  number  density,  m 
is  molecular  mass,  [O2]  is  the  oxygen  mole  frac¬ 
tion,  and  P„  and  T„  are  the  pressure  and  tempera¬ 
ture  of  the  ambient  gas. 

In  Fig.  2,  the  results  in  Fig.  1  have  been  cast 
into  the  form  of  collision  efficiency,  assuming  the 
reaction 


2C(s)  +  02^  2CO.  (5) 

Each  reactive  collision  with  O2  is  assumed  to 
remove  two  carbon  atoms;  therefore,  m  is  24  gm- 
atoms  for  this  reaction.  The  collision  efficiency 
monotonically  increases  with  decreasing  pressure. 
In  Fig.  3,  the  NSC  oxidation  rate  expression  is 


again  evaluated,  but  now  extrapolations  in  both 
temperature  and  pressure  are  shown.  The  lower 
oxygen  partial  pressures  in  Fig.  3  are  more  typical 
of  flame  conditions.  It  is  not  obvious  in  Fig.  3  that 
anything  is  amiss.  At  a  fixed  temperature,  the  oxi¬ 
dation  rate  is  a  monotonic  increasing  function  of 
oxygen  partial  pressure.  In  Fig.  4,  the  oxidation 
rates  in  Fig.  3  are  recast  as  collision  efficiencies. 
Over  this  pressure  range,  the  collision  efficiency  is 
not  a  monotonically  increasing  function  with 


1/Tempeiatute,  1/K 


Fig.  2.  Curve  fit  of  NSC  collision  efficiency. 


I/Temperature,  1/K 

Fig.  3.  Extrapolations  of  NSC  curve  fit  of  oxidation 
rates. 


Fig.  4.  Extrapolations  of  NSC  curve  fit  of  coiiision 
efficiency  versus  1/temperature. 
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decreasing  pressure.  As  a  convenience,  Fig.  5  pre¬ 
sents  the  same  collision  efficiencies  as  Fig.  4,  piot- 
ted  versus  temperature  rather  than  1 /temperature. 
Note  that  there  is  significant  overlap  of  the  curves 
in  the  1700  to  2200  K  temperature  range  -  an 
important  range  for  flames. 


Fig.  5.  Extrapolations  of  NSC  curve  fit  of  collision 
efficiency  versus  temperature. 


To  see  the  non-monotonic  behavior  more 
clearly.  Fig.  6  presents  the  collision  efficiency  as  a 
function  of  oxygen  partial  pressure  for  two  temper- 
atures-1250  and  2000  K.  For  1250  K,  the  collision 
efficiency  is  a  monotonically  increasing  function 
with  decreasing  pressure.  This  is  the  expected 
physical  behavior.  The  pressure  dependence  of  the 
oxidation  rate  is  related  to  the  occupation  of  active 
sites  on  the  soot  surface  by  other  oxygen  mole¬ 
cules  (or  atoms).  At  extremely  low  pressure,  the 
collision  efficiency  is  simply  the  ratio  of  the  active 
surface  area  to  the  total  surface  area.  This  ratio  is 
strictly  a  function  of  soot  temperature.  In  order  to  fit 
their  data,  Nagle  and  Strickland-Constable  used  a 
modification  of  the  surface  reaction  theory  of  Bly- 
holder,  Binford,  and  Eyring.®  This  theory  states  that 
the  number  of  active  sites  on  the  surface  is  a  func¬ 
tion  of  temperature  only,  and  is  therefore  indepen¬ 
dent  of  pressure.  Active  sites  are  occupied  by  oxy¬ 
gen  molecules  (or  atoms)  for  a  finite  time,  rendering 
these  sites  temporarily  inactive.  If  the  collision  fre¬ 
quency  is  high  enough,  the  molecules  occupying 
the  active  sites  will  interfere  with  other  molecules 
arriving  at  the  surface.  Therefore,  increasing  pres¬ 
sure  causes  an  apparent  decrease  in  the  number  of 
active  sites  on  the  soot  surface.  This  is  reflected  as 
a  decrease  in  the  collision  efficiency.  The  number 
of  active  sites  on  the  surface  is  fixed  by  the  temper¬ 
ature,  but  the  availability  of  an  active  site  is  affected 
by  the  collision  frequency.  At  high  collision  frequen¬ 


cies,  molecules  collide  with  the  soot  surface  faster 
than  the  occupied  sites  are  cleared.  This  means 
that  the  collision  efficiency  must  be  a  monotonically 
increasing  function  with  decreasing  pressure.  As 
seen  on  Fig.  6,  at  2000  K  the  collision  efficiency 
computed  from  the  NSC  expression  is  not  a  mono¬ 
tonically  increasing  function  with  decreasing  pres¬ 
sure.  The  efficiency  has  a  maximum  at  about  0.02 
atm  and  then  decreases  with  further  decrease  in 
pressure.  This  decrease  is  not  physical  and  vio¬ 
lates  the  Blyholder-Binford-Eyring  theory  on  which 
the  NSC  cunre  fit  is  based.  It  is  caused  simply  by 
extrapolating  the  NSC  expression  beyond  its  range 
of  validity. 


Fig.  6.  Pressure  dependence  of  NSC  collision 
efficiency. 


The  collision  efficiency  is  plotted  as  a  function  of 
pressure  for  various  temperatures  in  Fig.  7.  Note 
that  at  high  temperature,  the  collision  efficiency  is 
nearly  constant  with  pressure.  This  is  due  to  the 
extremely  short  residence  times  in  the  active  sites 
at  high  temperature.  The  pressure  at  which  the  col¬ 
lision  efficiency  is  maximum  is  the  lowest  valid 
pressure  for  extrapolating  the  NSC  rate  expression. 
The  locus  of  these  minimum  valid  pressures  is  plot- 


Oxygen  Partial  Pressure,  atm 

Fig.  7.  Pressure  dependence  of  collision  effi¬ 
ciency  for  various  temperatures. 
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ted  as  a  function  of  temperature  in  Fig.  8.  This  fig¬ 
ure  gives  a  rapid  estimate  of  the  pressure  at  which 
the  NSC  rate  expression  may  begin  to  be  suspect. 
However,  at  temperatures  greater  than  2500  K  the 
coilision  efficiency  is  relatively  flat,  so  the  error 
introduced  at  low  pressure  is  small. 


Fig.  8.  Minimum  valid  pressure  for  NSC  rate 
expression,  derived  from  maximum  colli¬ 
sion  efficiency. 

There  is  another  way  to  see  the  nonphysical 
behavior  of  the  NSC  rate  expression  without  con¬ 
verting  oxidation  rates  to  collision  efficiencies. 
Since  the  number  of  active  sites  on  the  surface  is  a 
function  of  temperature  only,  if  pressure  is  reduced 
by  some  factor,  the  oxidation  rate  must  reduce  by 
the  same  factor  or  less.  In  other  words,  if  the  pres¬ 
sure  is  reduced  a  factor  of  two,  the  oxidation  rate 
can,  at  most,  decrease  a  factor  of  two.  The  same 
obviously  holds  true  for  increases  in  pressure.  A 
doubling  of  pressure  can  produce,  at  most,  a  dou¬ 
bling  in  the  oxidation  rate.  Therefore,  a  simple 
check  is  to  normalize  the  pressure  and  the  oxida¬ 
tion  rate  by  their  values  at  some  given  pressure.  In 
Fig.  9,  the  normalized  oxidation  rate  is  plotted  as  a 
function  of  normalized  pressure  at  a  temperature  of 
2000  K.  The  pressure  is  normalized  by  0.5  atm, 


Fig.  9.  Normalized  oxidation  rate  versus  normal¬ 
ized  pressure  at  2000  K. 


and  the  oxidation  rate  is  normalized  by  the  oxida¬ 
tion  rate  at  0.5  atm.  At  a  normalized  pressure  of 
0.001 ,  the  normalized  oxidation  rate  is  0.0007.  This 
means  that  a  reduction  in  pressure  by  a  factor  of 
1000  produced  a  reduction  in  oxidation  rate  by  a 
factor  of  1500.  This  is  clearly  nonphysical  and  vio¬ 
lates  the  Blyholder-Binholder-Eyring  theory  on 
which  the  NSC  cunre  fit  is  based.. 

The  oxidation  rate  must  change  by  a  factor 
equal  to  or  less  than  the  factor  that  the  pressure  is 
changed.  If  the  oxidation  rate  changes  by  the  same 
factor  as  the  pressure,  then 

=  cons  tent  (6) 

Taking  the  logarithmic  derivative  yields 


or 


3<0 
<0  aPog 


=  1.0 


(8) 


Since  the  oxidation  rate  can  change  by  a 
smaller  factor,  the  actual  requirement  can  be 
defined  as 


Qs 


am 

"^aPog 


^1.0 


(9) 


The  indicated  partial  derivative  can  be  evalu¬ 
ated  analytically  from  Eq.  (1).  This  quantity  (G)  is 
plotted  as  a  function  of  pressure  for  various  temper¬ 
atures  in  Fig.  10.  Whenever  G  is  greater  than  one, 
the  oxidation  rate  is  changing  by  a  greater  factor 
than  the  pressure,  and  is  therefore  invalid.  Compar¬ 
ing  Figs.  8  and  1 0,  it  can  be  seen  that  G  rises  above 
unity  at  the  same  pressure  at  which  the  collision 
efficiency  is  maximum.  Therefore,  both  methods  of 
determining  the  minimum  valid  pressure  for  the 
NSC  rate  expression  at  a  given  temperature  give 
the  same  results.  The  locus  of  pressures  where  G 
just  rises  above  unity  is  plotted  as  a  function  of  tem¬ 
perature  in  Fig.  1 1 .  This  again  gives  an  indication  of 
the  minimum  pressure  for  which  the  NSC  rate 
expression  is  valid.  Note  the  similarity  with  Fig.  8. 
The  two  methods  for  determining  the  minimum 
valid  pressure  agree  almost  precisely. 
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Conclusions 


Fig.  10.  Pressure  dependence  of  G  for  various 
temperatures. 


Fig.  1 1 .  Minimum  valid  pressure  for  NSC  rate 
expression,  derived  from  G  >  1 . 

Impact  of  Extrapolation  on  Oxidation  Rates 

To  see  the  potential  impact  of  extrapolation, 
consider  a  flame  with  an  oxygen  partial  pressure  of 
0.001  atm  and  a  temperature  of  2000  K  (approxi¬ 
mating  conditions  of  Lin,  et  al.).  The  NSC  expres¬ 
sion  yields  a  collision  efficiency  of  8  x  1 0"^  at  these 
conditions.  However,  the  collision  efficiency  can  be 
no  less  than  that  at  the  lowest  pressure  of  the 
experimental  NSC  data  (0.10  atm)  which  is, about  3 
X 10'®.  Therefore,  the  NSC  expression  yields  a  col¬ 
lision  efficiency  that  is  at  least  a  factor  of  3-4  too 
low.  Of  course,  the  collision  efficiency  should  actu¬ 
ally  continue  to  increase  with  decreasing  pressure. 
If  the  collision  efficiency  at  2000  K  is  assumed  to 
have  the  same  pressure  dependence  as  the  effi¬ 
ciency  at  1250  K,  the  NSC  expression  at  0.001  atm 
may  be  low  by  as  much  as  a  factor  of  10  or  more. 
This  demonstrates  that  the  nonphysical  nature  of 
the  NSC  rate  expression  can  have  a  significant 
impact  on  soot  oxidation  rates  at  flame  conditions. 


The  NSC  rate  expression  is  invalid  when 
extrapolated  to  oxygen  partial  pressures  typical  of 
laboratory  flames.  The  extrapolations  exhibit  non¬ 
physical  behavior  and  produce  collision  efficiencies 
significantly  lower  than  would  be  expected.  The 
nonphysical  behavior  is  revealed  when  the  NSC 
oxidation  rate  expression  is  cast  in  the  form  of  colli¬ 
sion  efficiency.  Oxidation  rates  obtained  from  the 
NSC  oxidation  rate  expression  at  typical  flame  con¬ 
ditions  may  be  an  order  of  magnitude  lower  than 
expected  on  physical  grounds.  The  surface  reac¬ 
tion  rate  theory  of  Blyholder,  Binford,  and  Eyring 
should  not  produce  the  non-physical  behavior, 
since  the  theory  assumes  that  the  number  of  active 
sites  is  independent  of  pressure.  Future  work 
should  include  the  fitting  of  the  Nagle  and  Strick- 
land-Constable  data  using  the  unmodified  Bly¬ 
holder,  Binford,  and  Eyring  theory  to  ascertain  if 
the  modifications  introduced  by  Nagle  and  Strick¬ 
land  Constable  introduced  an  error.  A  new  rate 
expression  that  does  not  exhibit  this  nonphysical 
behavior  should  be  derived  and  compared  with  a 
variety  of  oxidation  rate  measurements  that 
encompass  a  wide  pressure  range. 
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pressure  is  reduced.  At  typical  flame  conditions,  extrapolation  of  the  NSC  curve  fits  yields  collision  efficiencies  that  may  be 
an  order  of  magnitude  too  low. 
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